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On  the  estimation  of  long  range  sound 
Intensities  by  ray  trace  methods 

C.  S.  Clay  . * 


The  problem  of  calculating  ray  trace  intensities  at  short  range  is 
relatively  simple  when  a few  particular  ray  paths  from  the  source  to  the 
receiver  are  known.  Our  interest  here  is  in  the  larger  ranges,  i.e.  , ranges 
much  greater  than  the  '^skip  distances^'.  Diffraction,  scattering  and  the 
changes  of  the  velocity  profiles  over  the  large  range  make  it  convenient, 
even  a necessity,  to  regard  the  ray  paths  as  being  distributed  over  a 
region  of  the  ocean. 

H t 

To  avoid  argument,  we  agree  that  the  particular  ray  paths  from  a 
source  position  to  a receiver  position,  regardless  of  the  ray,  can  be 
calculated  for  a sound  velocity  profile  over  the  path.  Temporal  changes 
and  lack  of  exact  knowledge  of  the  sound  velocity  profile  make  this 
deterministic  problem  a futile  exercise.  uae  the  ray  tracing 

W.n  t - vr.'.V 

procedure* primarily  to  determine  the  aperture  at  the  source  and  the 
distribution  of  the  energy  in  the  region  of  the  receiver.  Hence,  we  eau 
only  vEflmifitv  the  average  sound  intensity,*"  '"  e>‘  ' • |v 

On  Fig.  1,  wo  have  shown  several  different  classes  of  rays.  At 
short  range,  one  has  the  familiar  ray  diagrams.  At  large  range,  the  dis- 
tribution of  the  rays  produces  a.Jtoid  somewhat  like  that  indicated  on  the 
right.  The  angular  aperatures  of  irisoniflcation  are  AG  A 8^  and  AG(  . 

For  a surface  duct,  the  attenuation  limits  the  range  to  short  distances  and 
most  of  its  energy  eventually  goes  into  leakage  arrivals.  Their  contribution 
Is  generally  small.  At  short  range,  the  hydrophones  "a"  and  "b"  may  or 
may  not  observe  signal  as  dependent  upon  the  insonification.  At  large  range, 
the  energy  Is  distributed  between  the  upper  and  lower  reflection  or  turning 
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Fig.  1 Types  of  ray  paths.  The  areas  of  insonification  are  shaded  and 


onlv  down  aoina  ravs  are  shown. 


depths,  "c"  and  "d"  would  receive  the  signal  for  all  cases.  When  "g"  is 
beneath  the  turning  depth,  as  shown  for  the  refraction  case,  it  would  not 
receive  signals  except  by  scattering  and  leakage  effects. 

In  the  calculation  of  the  intensities  to  follow,  we  assume  that  the 
intensities  due  to  each  type  of  propagation  can  be  determined  separately 
and  then  added  together.  The  diagram  on  Tig.  2 illustrates  the  geometry  for 
the  refracted  arrivals. 


Tig.  2 Maximum  and  minimum  rays  for  the  source 
and  receiver  depths  and  z^ 

From  the  figure,  it  is  evident  that  the  energy  between  the  rays  from  the 

source  is  distributed  over  the  depth  at  the  right.  Let  us  assume  that  the 

point  source  gives  an  rms  sound  pressure  p^  at  distance  r^  from  the 

source.  Ignoring  the  negligible  change  of  c(z)  near  the  source  and  for 

small  A0,  the  outgoing  energy  is  proportional  to  2 tt  pQ  r^  cos  0 

(0  - 0 ).  Ignoring,  for  the  present,  the  absorption  and  other  losses, 

max  min  _ 2 

the  energy  at  distance  r is  proportional  to  2 rr  r \ p dz. 
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lire  constant  of  proportionality  is;  the  relation  between  intensity  and 
\'  p“  dz  and  Is  nearly  the  same  tor  both  the  source  and  receiver.  By  con- 
servation ot  energy , we  can  write 


ll  n r h 2 

p r cos  0 6 r \ p dz 

O O p 
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(l) 


se  * e - 0 

max  min 

0 - (0  * 0 )/2 
max  min 

h » depth  of  the  ocean  (~) 

where  „\0  is  in  radians  and  i , r,  and  h aie  m the  same  units.  Sometimes 

o 

the  lange  at  which  the  spieadiiui  changes  from  spherical  to  cylindrical 

spieadlng  is  used  to  define  a quantity  equivalent  to  .\0  (Urick,  h)t>7). 

) 

Hie  dependence  of  p“  upon  depth  is  an  important  quantity  in  the 

detei  mi  nation  of  the  signal  at  a receiver.  At  the  beginning,  we  suggested 

that  it  th«'  receiver  is  above  the  minimum  turning  depth  or  beneath  the 

maximum  turning  depth,  the  field  is  zero.  this  kind  of  thinking  is  helpful 

in  making  qualitative  descriptions  ot  the  problem,  however  some  better 

means  is  needed  to  estimate  p-  as  a function  of  depth.  As  has  been  done 

in  the  past,  we  suggest  letting  the  intensity  ot  rays  passing  through  a 

> 

depth  be  proportional  to  p“  . 

Motor e we  begin,  let  us  make  a tew  comments  concerning  tin' 
wave-ray  ot  paiticle  dualism  ot  physics.  Hectrons  are  sometimes  treated 
as  waves  and  sometimes  as  particles,  l‘he  formalism  is  chosen  to  be 
convenient  to  the  problem.  Hie  same  sort  of  dualism  applies  in  classical 
physics  to  light  and  sound  waves.  It  the  anertures  are  very  large  and  the 
wave  fronts  are  nearly  plane  over  many  wave  lengths,  then  it  is  both 


convenient  and  accurate  to  use  the  rays  or  normals  to  the  wave  fronts  to 
calculate  travel  times  and  intensities.  Wave  descriptions  are  needed  when 
the  aperture  is  a small  number  of  wave  lengths.  Also,  wave  descriptions 
are  needed  whenever  rays  are  focused  because  one  is  looking  at  an  inter- 
ference phenomenon.  Generally  in  these  situations,  one  should  say  that 
a wave  field  is  in  the  region  of  the  receiver  rather  than  that  a ray  or  bundle 
of  rays  has  passed  through  the  receiver. 

At  large  range,  the  wave  front  associated  with  a ray  has  spread 

enormously  and  thus  has  had  many  interactions  with  the  upper  and  lower 

boundaries  of  the  ocean.  The  surface  and  bottom  limit  the  aperture  of  the 

transmission  system  and  cause  the  energy  to  be  spread  over  a range  in 

depth.  We  will  not  attempt  to  calculate  the  spreading  of  the  sound  energy 

but  will  assert  that  it  is  reasonable  to  use  the  average  ray  density  to 
2 1 

estimate  p . 

As  shown  on  Fig.  3,  the  ocean  at  range  r has  been  subdivided  into 

small  regions  bounded  by  r,  . . . r . . . , r. , and  z , z,  , . . . z.T  . The 

1 m Mol  N 

ray  traces  from  the  source  are  made  at  equal  angle  increments.  A number  of 

rays  have  been  shown  as  passing  through  the  regions.  In  keeping  with  our 

view  that  the  wave  field  associated  with  the  ray  is  rather  large,  we  want 

to  average  the  ray  intensity  over  the  distance  Ar  for  the  depth  increments 

z , z,  , etc.  By  ray  tracing,  the  number  of  rays  passing  through  the  area 
o 1 

at  r , (z  _ - z ) is  P . The  average  number  of  rays  passing  through 
m n+ 1 n m , n 

the  depth  increment,  z , - z , is 
K n+1  n 

M 

P e -J-  y p (3) 

n M L m,n 
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^At  the  shorter  ranges  it  is  desirable  to  include  the  x-z  spreading  factor. 
M.  Holl  (1967). 


Fig.  3 Rays  passing  through  small  regions:  The  area  at  r , (zn+j  “ zn) 
is  darkened. 


The  total  number  of  rays  N is 


N 

Total  number  of  rays,  P 5 Y Pn 

o 


(4) 


(Note  that,  if  no  rays  are  lost  in  Ar  , then  P is  the  sum  of  the  rays 
passing  through  the  wave  duct.) 

On  the  basis  of  our  assumption  that  the  mean  square  acoustical 

pressure  p2  at  the  depth  (z  + z 0/2  is  proportional  to  P thus  \ p2  dz 
K — n n n+ 1 n j 

and  p^  are  as  follows: 
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p“  ^ p ( cos  0-6  r “ r 1 (z  - z ) 1 (6) 

n o \ P / o n+ 1 n 


In  the  computation  of  P , the  averaging  distances  should  be  chosen  on  the 

basis  of  the  expected  detail  of  the  sound  field  and  the  range.  At  very 

large  range,  it  is  suggested  that  Ar  be  at  least  a "skip  distance"  or  more 

than  60  km  and  that  the  (z  - z ) be  many  sound  wave  lengths.  At  short 

n+ 1 n 

range,  Ar  should  be  small  enough  to  reveal  convergence  zones  and  shadow 

zones.  Specification  of  analytical  or  quantitative  values  of  Ar  and 

(z  - z ) is  difficult  and  beyond  the  scope  of  this  note. 
n+ 1 n ^ 

Equation  (6)  can  be  used  for  any  type  of  propagation.  For  example, 

the  amount  of  energy  trapped  in  the  surface  duct  is  given  by  and  its 

corresponding  insonified  depths.  The  effective  sources  for  leakage  signals 

can  be  approximated  as  being  at  the  bottom  of  the  duct  and  it  is  reasonable 

to  assume  the  ray  directions  are  slightly  downward. 

V 

Deep  water  ray  trace  solutions  often  have  rays  that  have  had  a few 
bottom  interactions  along  the  path.  A way  of  treating  these  rays  is 
suggested  on  Fig.  lb.  Some  of  the  energy  is  lost  at  each  bottom 
reflection  and  part  is  reflected.  At  a reasonably  large  distance,  the 
energy  is  again  distributed  between  the  maximum  and  minimum  depths. 

Thus,  we  can  include  n bottom  reflections  as  follows: 

y pbr2  - % c°s  v aeb,  ro2  w~‘  ri2  r22  • • • r„2  <7) 

where  R , , . . . Rn  are  the  average  reflection  coefficients  along  the 

I : 
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transmission  path . Obviously,  if  the  bottom  is  deep  enough  over  part  of 
the  path  such  that  the  sound  turns  above  the  bottom,  one  only  need  include 
the  losses  where  the  sound  reflects  at  the  bottom.  The  total  sound 
intensity  is  the  sum  of  the  contributions  of  refracted,  surface  and  bottom 
reflected  paths. 

Although  it  is  obviously  a poor  practice,  sometimes  one  must  assume 
a density  distribution  of  the  rays.  If  one  assumes  that  the  distribution  is 
linear  from  the  bottom  to  top  turning  points,  then  Eq.  (6)  can  be  written  as 
follows: 


pZ  (z  , - z ) P“  (z  - z . )/l 

n ml  n max  nun 


r N = P 
n 


— r 2 Z - 1 - 1 

p~  p “ cos  6 .X  6 r “ (z  - z ) r 

o o max  nun 


But,  we  do  not  suggest  the  use  of  Eg.  (10)  if  rough  estimates  of  the  density 
are  available. 

Tor  the  convenience  of  the  reader,  the  transmission  loss  in  dB  is 
with  the  aid  of  Eq.  (t>)  and  including  sound  absorption,  the  following: 


TL  * 10  log  r + Ar  - 10  log  (P  /P) 

n n 


+ 10  log  (z  - z ) - 10  log  cos  6 

n+1  n 


10  log  A 6 - 10  log  r “,  dB 

o 


where 


A is  the  sound  absorption  in  dB/m 


P is  the  average  number  of  rays  passing  through  the  depth 
n 

increment  (z  j - z^) 

P is  the  total  number  of  rays  at  distance  r 

A9  is  the  angular  aperture  of  the  source  in  radians 

r z and  r are  in  the  same  units,  m. 

' n'  o 
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